THE SPECTRAL DENSITY OF THE SCATTERING MATRIX FOR 

HIGH ENERGIES 
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Abstract. We determine the density of eigenvalues of the scattering matrix of the 
Schrodinger operator with a short range potential in the high energy asymptotic regime. 
We give an explicit formula for this density in terms of the X-ray transform of the potential. 



1. Introduction and Main Result 

1.1. Introduction. The object of study in this paper is the (on-shell) scattering matrix 
S(X) corresponding to the scattering of a d- dimensional quantum particle on an external 
short range potential V at the energy A > 0. The scattering matrix 5(A) is a unitary 
operator in L 2 (E> d ~ l ) and the difference S(X) — I is compact. Thus, the eigenvalues of S(X) 
can be written as 

(1.1) exp(i0„(A)), 9 n (X) e l-n, n), neN 

and & n (X) — > as n — > oo. The quantities 9 n (X) are known as scattering phases. The scat- 
tering phases are usually discussed in physics literature (see e.g. [BJ Section 123]) under the 
additional assumption of the spherical symmetry of V; here we do not need this assumption. 

Our aim is to study the asymptotic distribution of scattering phases {6 l n (A)}^ =1 when 
A — > oo. It turns out that after an appropriate scaling the asymptotic density of scattering 
phases can be described by a simple explicit formula involving the X-ray transform (see 
( II. 5p ) of the potential V. This formula has a semiclassical nature. 

The key idea of this paper goes back to the work of M. Sh. Birman and D. R. Yafaev [T] 
(see also [2]), where the asymptotics of 9 n (X) for a fixed X and n — > oo was determined for 
some class of potentials V. This asymptotic behaviour is not uniform in A, and thus our 
results cannot be derived from those of pp. However, both the results of [1] and our results 
are based on the following observation. The leading term of the asymptotics of 6 n (X) (in 
both asymptotic regimes) is determined by the Born approximation of the scattering matrix. 
The Born approximation is essentially a pseudodifferential operator (\l/DO) on the sphere 
E> d ~ l with the symbol given by the X-ray transform of V. Standard \l/DO results can be 
used to give spectral asymptotics for operators of such type. In both asymptotic regimes, 
the desired spectral asymptotics are given by a Weyl type formula involving the symbol of 
the tfDO. 

1.2. The scattering matrix. Let us briefly recall the relevant definitions. Let Hq = —A 
and H = — A + V be the Schrodinger operators in L 2 (M. d ), d > 2, where V is the operator of 
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multiplication by a real- valued potential V G C(lR d ) which is assumed to satisfy the estimate 
(1.2) <(T^. p>i 

with some constant C > 0. It is one of the fundamental facts of scattering theory [U [5] that 
under the assumption (jl.2p the wave operators 

W± = s-lim e itH e^ tHo 

t— >±oo 

exist and are complete; the scattering operator S = W+W- is unitary in L 2 (IR d ) and com- 
mutes with H . Let F : L 2 (M d ) -> L 2 ((0, oo); L 2 (S d-1 )) be the unitary operator 

(Fu)(X,u) = ^\( d ~^ /4 u(V\u), A>0, ueS d -\ 
v2 

where u is the usual (unitary) Fourier transform of u. The operator F diagonalises Ho, i.e. 

(FH u)(X,uj) = X(Fu)(X,u), Vu G C^(R d ). 

The commutation relation Sifo = HqS implies that F also diagonalises S, i.e. 

(FSu)(X,-) = S(X)(Fu)(X,-), 

where S(X) : L 2 (E> d ~ 1 ) — > L 2 (S a!_1 ) is the unitary operator known as the (on-shell) scattering 
matrix. See e.g. the book [9] for the details. 

Under the assumption (11.21) the operator S(X) — I is compact and consequently the eigen- 
values of S(X) (enumerated with multiplicities taken into account) can be written as (11.11) 
with 9 n (X) — > as n — > oo. 



1.3. The purpose of the paper. Our purpose is to describe the asymptotic density of the 
eigenvalues of 5(A) as A — > oo. We recall the estimate (see e.g. [9., Section 8.1]) 

(1.3) ||5(A)-/|| =0(A~ 1/2 ), A^oo. 

This estimate is sharp, i.e. 0(A -1 / 2 ) cannot be replaced by o(A -1 / 2 ); this can be seen 
by considering the case of a spherically symmetric potential and using the separation of 
variables. Thus, the spectrum of the scattering matrix S*(A) for large A consists of a cluster 
of eigenvalues located on an arc of length 0(A -1 / 2 ) around 1. Let us define the eigenvalue 
counting measure for S(X) . The estimate (jl.3p suggests the following scaling: for A > 1, set 
k = a/X > and define (using notation (II. ip ) 

(1.4) MM) = #{n e N : a < k6 n (k 2 ) < /?}, [a, 0\ C M\{0} 

where # denotes the number of elements in the set. We will study the weak asymptotics of 
Hk as k — > oo, i.e. we consider the asymptotics of the integrals 




i/j(t)dfj,k(t), k — > oo 



for test functions ip G C °°(M\{0}). 
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1.4. Main result. In order to describe the weak limit of the measures we need to fix 
some notation. For any u G § d_1 , let A w C M. d denote the hyperplane passing through the 
origin and orthogonal to u. We equip both and with the standard (d— l)-dimensional 
Lebesgue measure (=Euclidean area). We set 

1 f°° 

(1.5) X{u, r}) = -- V(tu + rj)dt, to E S d ~\ V e A„. 

^ J -oo 

The function X (up to a multiplicative factor) is known as the X-ray transform of V in 
the inverse problem literature. The following elementary estimate is a direct consequence of 

(1.6) \X(",V)\ <C(V0(1 + M) 1_P , ue§ d -\ t|G A u 
with some constant C(V). We define a measure \x on M\{0} by 

fi([a,0]) = (27r) 1 - d meas{(a;,r/) e S^ 1 x A w : a < X{uj,rj) < /?}, [a, /3] C M\{0}, 

where meas denotes the usual product measure. By the boundedness of V, the measure \i 
has a compact support. The measure \i need not be absolutely continuous. The measure 
/i may be weakly singular at zero in the following sense: p((0, oo)) or oo,0)) may be 
infinite, but, by the estimate (11 .6p we have 

/oo 
\tfdn(t) < oo, W>(d-l)/(p-l). 
-oo 

Our main result is as follows: 

Theorem 1.1. Let V G C(lR d ) be a potential satisfying ( 11.2H . TTien /or any test function 
feC »(R\{0}), 

/OO /* oo 

V»(*)d/x fc (0 = / ^{t)d/i{t). 
-oo J-oo 

This can be more succinctly put as the weak convergence of the measures 

(1.9) k x ~ d \i k -» p, -)■ oo. 

Much of the inspiration for both the content of this paper and the proofs may be found in 
[7], where similar asymptotics are determined for the spectrum of the Landau Hamiltonian 
(i.e. two-dimensional Schrodinger operator with a constant homogeneous magnetic field) 
perturbed by a potential which obeys the same condition (jl.2p . 

We would like to mention also the paper [TO] where the high energy asymptotic distribution 
of the phases 6 n (X) was studied for scattering problems on manifolds of a certain special class. 
The results of [10] are much more detailed than ours and include the asymptotics of the pair 
correlation measure. 

1.5. Comparison with pp. In [TJ, the case of potentials with the power asymptotics at 
infinity of the type 

(1.10) V(x) = v(x/\x\)\x\- p (l + o(l)), \x\-Kx, p>l, 
was considered. Using our notation p^, p, the result of [1] can be written as 

^-^((e, oo)) ~ //((e, oo)), 

k l ~ d jJL k ((-00, -£)) ~ //((-CO, -£)), 
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when k > is fixed and e — > +0. Here a ~ b means | — > 1. 

Clearly, our main result (11. 9ft is expressed by the same formula as ( ll.lip . but the asymp- 
totic regimes are different. Neither of the results (jl.9p . ( ll.lip implies the other one. 

1.6. Semiclassical interpretation. By the definition of the scattering operator S, for any 
if) G L 2 (R d ) we have 

i((S - I)ifj,if)) = » Hm((e- 2i * H e^ V,e- i * Ho V) - ||Vf ) = * / -f (e-^e^, e- itH °V)dt 

/•oo 

(Ve" 2 ^'^, e- itHo ip)dt + / (Ve** ^, e^e^^dt 

Jo 

If ^ corresponds to large energies, the right hand side can be approximated by the first term 
in its expansion in powers of V. This means that we can replace e ttH by e ttH ° in the above 
expressions, and so 

/oo 
(Ve- itHo ip, e- itHo ip)dt, V e ^ 2 (K d ), 
-oo 

which is exactly the Born approximation in the time- dependent picture. 

In order to write down the classical analogue of the right hand side of (I1.12p . assume 
that if) is concentrated near x in the coordinate representation and near p in the momentum 
representation. Then if) represents a particle with the coordinate x and momentum p, and 
in the same way e~ ltH °ip represents a particle with the coordinate x + 2pt and momentum 
p. Thus, the classical analogue of the right hand side of ( I1.12p is 



/OO /*oo 
V(x + 2pt)dt = — / V(x + ut')dt', (x,p) E 
2|P| J-oo 



where u = & G E> d 1 . This calculation explains the appearance of the X-ray transform in 
the asymptotics of S(X). 

1.7. Key steps of the proof. First we recall the stationary representation for the scattering 
matrix. For k > and p > 1, we define the operator T p (k) : L 2 (R d ) — > L 2 (S d_1 ) by 

(1.13) (T p (k)u)(u) = -^k {d ~ 2)/2 (2n)- d/2 [ u{x){x)- p/2 e- lk{x ' uj) dx 1 

V 2 JE d 

where (x) = (1 + l^ 2 ) 1 / 2 . By the Sobolev trace theorem, T p (k) is bounded for p > 1. Next, 
let 

T{z) = (x)- p/2 (H - zl)- l {x)- p/2 , lm 2 >0; 

according to the limiting absorption principle, the limits T(k 2 ± iO) exist in the operator 
norm for all k > 0. Denote by J the bounded operator of multiplication by (x) p V(x) 
in L 2 (M. d ). The stationary representation for the scattering matrix can be written as (see 
e.g. Section 6.6]) 

S(k 2 ) = I - 2mT p (k)(J - JT(k 2 + iO)J)T p (k)*, k > 0. 

The asymptotic density of eigenvalues of S(k 2 ) for large k is determined by the Born ap- 
proximation of the scattering matrix, defined as 

(1.14) S B (k 2 ) = I -2mT p (k)JT p (k)\ k > 0. 
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The key observation due to M. Birman and D. Yafaev [J] is that the operator T p (k) JT p (k)* 
in L 2 (S d_1 ) with the integral kernel 

JR d 

can be represented as a \l/DO on the sphere with the symbol (up to inessential constants) 
X(u,rj). We combine this observation with the standard semiclassical pseudodifferential 
techniques to obtain the spectral asymptotics of Ssik 2 ). In this way we prove the asymptotic 
formula (see Lemma [3.21) 



/oo 
t e dp(t) 
-oo 



for any natural number £ > (d — 1) / (p — 1); note that the r.h.s. of fl 1 . 1 5 [) is finite by (11.71) . 

Using (11.151) and the estimates for the Schatten norm of S(k 2 ) — Ssik 2 ) we prove that 
(11. 8p holds true for test functions if)(t) which coincide with t e for all sufficiently small t. 
Theorem 11.11 follows by an application of the Weierstrass approximation theorem. 

1.8. Acknowledgements. The authors are grateful to Y. Safarov for expert advice con- 
cerning the \l/DO aspect of this work. The authors are grateful to N. Filonov and D. Yafaev 
for a number of useful discussions and remarks on the text of the paper. 



2. Preliminary statements 

2.1. The limiting absorption principle and its consequences. First we need some 
notation. We denote by S p , p > 1, the usual Schatten class and by || • \\ p the associated 
norm. Let X p be the normed linear space of all potentials V G C(R d ) satisfying (TO} with 
the norm 

\\V\\ Xp = su V \V(x)\(xy. 
We recall the following estimates from [9]: 

Proposition 2.1. Let V G X p with some p > 1. Then for any £ > 1 satisfying £ > one 
has 

(2.1) snpk^WklmSBik 2 )^ < C(£, p, d) \\V\\ Xp , 

k>l 

(2.2) supA^H^Im^A; 2 ) - S(k 2 ))\\t < C(£,p,d,V), 

k>l 

(2.3) supA^HATm^A; 2 )^ < C{l,p,d,V). 

k>l 

The estimate ( 12.11) is a direct consequence of (HJ Proposition 8.1.3]. The estimate ( 12. 2h 
is proven in P, Proposition 8.1.4]. The estimate ( 12. 3 h is a direct consequence of ( 12.11) and 

Lemma 2.2. Let V G X p with p > 1. Then for any integer £ > 1 satisfying £ > dj ^, one 
has 

| Tr(klmS{k 2 )Y - Ti(klmS B (k 2 )Y\ = 0(k d - 2 ), k oo. 
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Proof. From 

e-i 

A e -B e = J2 A K A ~ B)B e - 1 -J 

3=0 

one easily obtains 

(2.4) \Tr(A e -B*)\<e\\A-B\\ e max{\\A\\*-\\\B\\ e f 1 }, A,BeS e . 
Thus, it suffices to prove the relation 

\\klm{S{k 2 ) - S B {k 2 ))\\ e max{\\kImS{k 2 )\\ e e -\ || jfclm S B {k 2 ) \\ e f 1 } = 0{k d ~ 2 ), k oo. 
The latter relation follows by combining (I2.ip - fl2.3p . □ 

2.2. Semiclassical $DO on the sphere. A semiclassical \l/DO in L 2 (S d_1 ) can be defined 
in a variety of ways; below we describe a slightly non-standard approach to this definition, 
which will simplify our exposition in Section |3j 
For u, u' G S d_1 such that u + lo' ^ 0, we set 

(2.5) k = k(u, lo') = U + U \ G S d -\ 

\CO + Co'\ 

Clearly, k(u,lo') is a smooth function of (u, lo') G x S^" 1 away from the anti-diagonal 

AD = {(lo, u') G S^ 1 x S^ 1 | u + lo' = 0}. 

In order to overcome the (inessential) difficulties related to the singularity of k at the anti- 
diagonal, we will assume that our amplitudes vanish in an open neighbourhood of AD. We 
will say that a function b = b(co,Lo',r)), to, to' G S d_1 , rj G A^^n, is an admissible amplitude, 
if: 

• b is a C°°-smooth function of its arguments; 

• b(u,u',r)) = if \t)\ is sufficiently large; 

• b(u,u',i]) = if (to,u') are in an open neighbourhood of AD. 

For an admissible amplitude b and a semiclassical parameter k > 0, we define the operator 
Op fc [6] in L 2 (S d_1 ) via its integral kernel 

(2.6) Op Jk p»](w,w / )=fAy / e -^-^>fc(u;,u/,77)d77, 

where w, a;' G § d_1 . It is easy to see that for lo ^ to' one has 

Op k [b)(u,u r ) = 0(k-°°), k^oo. 

This shows that the values of the amplitude b(no, to' , 77) away from an open neighbourhood 
of the diagonal lo = lo' do not affect the asymptotic properties of the operator Op fc [6] as 
k — > 00. 

Proposition 2.3. For any admissible amplitude b and any k > 0, the operator Op k [b] is 
trace class, and for any leN one has 

limf^-) + Ti(Op k [b]Y= [ [ b(u,u, V Yd V du. 
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Sketch of proof . The proof follows standard methods of \l/DO theory; see e.g. Theo- 
rem 9.6] for a similar statement in the context of operators in M ra . Here we only outline the 
main steps. 

First note that for £ = 1 the result of Proposition 12.31 is trivial, since by a direct evaluation 
of trace we have the identity 

k ^ ~ d+1 



(2.7) ^—j Tr(0 Pfc [6]) = ^ b(cu,u, V )d V du. 

Next, using the local coordinates on the sphere and the composition formula for symbols of 
^DOs (see e.g. [3] Proposition 7.7]), we obtain the following statement. For any N > 
there exists M > such that (Op k [b]Y can be represented as 

M 



(2.8) (Op k [b]Y = J2 k ~ m °^ b ^ + Rh 



m=0 



where b m are admissible symbols, bo is such that 
(2.9) b (u,u,r]) = b(u,u,r)Y, Vw G S^ 1 , V77 G A u , 

and Rk is an integral operator with a smooth kernel which satisfies 

sup|i?fc(o;, oj')\ = 0(k~ N ), k — > 00. 

LU,LU' 

Now taking N > d — 1 and evaluating the trace in (12. 8p . we get 

lim ( f) Tr(0 Pfe [6])^ = lim Tr(Op fc N) = / / b (u,u,r ] )dr ] dco. 

In view of (I2.9p . this proves the required identity. □ 



3. Proof of Theorem 11.11 
3.1. The Born approximation with V G C™(R d ). 
Lemma 3.1. Let V G C^(R d ). Then for any £ G N, one has 



(3.1) lim k 1 ~ d Ti(k\m.S B {k 2 )Y = I t e dfi(t) 



Proof. 1. For ease of notation we write Q{k) = klm Ssik 2 ). By HI. 131) and (11.141) . Q(k) is 
the integral operator in L 2 (S d_1 ) with the integral kernel 



Q(k)(co,u') = -U^-) [ e-^^'^V(x)dx, cu,cj'eS d - 1 . 
2 \2irJ J Rd 

For fixed u, u', u> + u' ^ 0, let k = k(u, uj') be as in H 2 . 5 f) . Write any x G M d as x = nt + 77 
with t G M, 7] G A K . Note that by the orthogonality relation (a; — u/) _L k, one has 

(u) — u', x) = (oj — u/, rf). 
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Thus, the integral kernel of Q(k) can be rewritten as 



If k^- 1 



(3.2) Q(k)(u,<J) = --l — ) I I e-^—'^ViKiu^y + ^dtdv 



oc 



2Ti 



e -*(»-<s*)X(K{u,d),Tj)d.Ti, 

where X is given by (II. 5p . From here we directly obtain the required identity ( 13 .ip in the 
case £ = 1 by integrating the kernel of Q(k) over the diagonal. Now it remains to prove (13. ip 
for £ > 2. 

2. Let xo £ C^°(E> d ~ l x § d_1 ) be such that xoiWjU)') = 1 in an open neighbourhood of the 
diagonal uj = uj' and Xo(^, w') = in an open neighbourhood of the anti-diagonal w + u/ = 0. 
Denote xi — 1 — Xo, and let 

Q{k) = Q {k) + Q 1 {k), 

where Qj{k) is the operator with the integral kernel Xj{w,w')Q(k)(u,u'). By the fast decay 
of the Fourier transform of V and by the fact that \u — uj'\ is separated away from zero on 
the support of xi, we see that 

swp\Qx(k)(uJ,L)')\ = 0(k~°°) } k^oo. 

Ul.UJ 1 

Thus, it suffices to prove that 

/CO 
t e dp(t). 
-00 

From (13.21) it follows that Qo(k) can be represented as a semiclassical \l/DO on the sphere of 
the type (1231) : 

b is an admissible amplitude in the sense discussed in Section [2T2| and X is given by (11. 5p . 
Applying Proposition 12.31 we get 

(3.3) lim fc 1 - d Tr(Q (fc))^ = (27r) 1 - d / / X(uj,r]Ydr] duj = [ t £ d/j(t), 

k ^°° Js d -i Ja„ V-co 

as required. □ 
3.2. The Born approximation with general V. 

Lemma 3.2. Let V G X p with p > 1. Then for any integer £ > 1 satisfying £ > 2=1 one /ias 

/oo 
t £ d/x(t). 
-oo 

Proq/. Let X° be the closure of C£°(R d ) in X p . For any £ > ^1, denote 

/oo 
t*d M (t), 
-oo 

gj{V) = limsupA^Tr^ImS^A; 2 )^, 



fc— >-oo 



</7(y) = liminf fc^TrtfclmSi^fc 2 ))' . 

fc— >co 
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By Lemma EU for all V G C£°(M d ) we have 
(3-4) 9e(V) = gf(V) = 9e{V). 

Recall that >Sb(£; 2 ) depends linearly on V . Using this fact and the estimates (12. ip and (12. 4p . 
it is easy to check that gf(V) are continuous functionals on X p . Next, using the last equality 
in (I3.3P and the estimate d 1 . 6 1) . it is easy to see that ge(V) is a continuous functional on X p . 
Thus, by an approximation argument, H3.4[) holds for any V G X p . Finally, for any V G X p 
and a given £ > ^5j> choose p\ such that 1 < p\ < p with £ > -^zj- Then X p C X pi and the 
previous argument proves (13.41) for all V G X pi which suffices. □ 

3.3. From the Born approximation to the full scattering matrix. 

Lemma 3.3. Let V G X p with p > 1. Then for any integer £ > 1 satisfying £ + 2 > 



lim k^ d / t e dp k {t) = 

k-^oo I 

J — oo 

Proof. By Lemmas 12.21 and 13. 2\ it suffices to prove 



t e dp(t). 



(3.5) 



lim k x - d 

k—^oo 



t e dp k (t) -Ti(lmkS(k 2 )Y 



0. 



Recalling the definition (11.41) of the measure p k , one sees that (13.51) is equivalent to 



(3.6) 



lim k 1 -^ 

fc— >oo 



J2l( 6 n(k 2 )Y ~ (sin 9 n (k 2 )Y 



n=l 



0. 



By (II. 3p we have < \0 n (k 2 )\ < tt/4 for all sufficiently large k and all n. From the elementary 
estimates \9 n \ < 2| sin ^ | and \9 n — sin# n | < C\ sin^] 3 which hold for \9 n \ < 7r/4, it follows 
that for all sufficiently large k 

oo 

k i- d+ iJ2\(e n y - (sm0 B )'| < k l - d +"Y, ( \°n ~ sm^nl £ \0 n \*\ sh^- 1 "'" 

oo 

< ^-^C^^lsin^, 



n=l 



n=l 



j=0 

^(£)||Im^(F)||; +2 . 



n=l 



Now (13. 6p follows by applying the estimate (I2.3P for ||Im/cS'(/c 2 )||£ + 2 to the result just ob- 
tained. □ 

Proof of Theorem \l.l[ By the estimate (ll.3p . the supports of p k are bounded uniformly in 
k > 1. On the other hand, by the boundedness of 1/, the support of p is also bounded. 
Thus, we may choose T > such that 

supp/i C [— T, T] and supp p k C [— T, T] for all > 1. 



Next, fix ip G C^°(1R\{0}), and let £o be an even natural number satisfying £ > — j-. Since 
ip(t) vanishes near t = by assumption, the function ip(t)/t io is smooth. By the Weierstrass 
approximation theorem, for any e > there exists a polynomial V'o(^) such that 

m)t- e °-Mt)\<z, VtG[-T,T]. 
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Denoting ip±(t) = (ipo(t) ± e)t e °, we obtain 

(3.7) V-W < m < Mt), Vt e l-T,T], 

(3.8) ip+(t) -ip-(t) = 2et e ". 
By (J3TD, we get 

POO POO POO 

(3.9) / ^{t)dfi k {t) < / i(){t)dfi k {t) < / 7P + (t)dfx k (t). 



By construction, ip±(t) are polynomials which involve powers t m with m > £ . Thus, we can 
apply Lemma [3731 to (13. 9p . which yields 



lim sup A; 1 ^ / ip(t)dfx k (t) < / V+(*)dju(*)> 

fc— >oo J —oo J —oo 



/oo /»oo 
ip(t)df* k (t) > / V-(*)d//(t). 
oo J-oo 

On the other hand, by (1377D . (TO . 

/*oo /*oo 

V-(t)d//(t) < / ^(^)d/i(t) < / V+(*)d/x(*) 



and 

/oo /*oo POO 

i/j+(t)dn(t) - / ip-(t)dn(t) = 2e / t £ °d/i(t). 
-oo J —oo J —oo 

By ( 11. 6ft . the integral in the right hand side of the last estimate is finite; denote this integral 
by C. Combining the above estimates, we obtain 

/oo poo 
ip(t)dfM k (t) < / il>(t)dn(t) + 2eC, 
-oo J —oo 

/OO POO 
Mt)dfi k (t) > / ip(t)dfi(t) - 2eC. 
oo J-oo 

Since e > can be taken arbitrary small, we obtain the required statement. □ 
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